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Abstract 

In statistical physics, the XY model in two dimensions pro- 
vides the paradigmatic example of phase transitions medi- 
ated by topological defects (vortices). Over the years, a 
variety of analytical and numerical methods have been de- 
ployed in an attempt to fully understand the nature of its 
transition, which is of the Berezinskii-Kosterlitz-Thouless 
type. These met with only limited success until it was real- 
ized that subtle effects (logarithmic corrections) that mod- 
ify leading behaviour must be taken into account. This re- 
alization prompted renewed activity in the field and signifi- 
cant progress has been made. This paper contains a review 
of the importance of such subtleties, the role played by vor- 
tices and of recent and current research in this area. Direc- 
tions for desirable future research endeavours are outlined. 



1. PHASE TRANSITIONS 

Phase transitions are amongst the most remarkable and ubiq- 
uitous phenomena in nature. They involve sudden changes 
in measurable macroscopic properties of systems and are 
brought about by varying external parameters such as tem- 
perature or pressure. Familiar examples include the tran- 
sitions from ice to water, water to steam, and the demag- 
netization of certain metals at high temperature. These 
dramatic phenomena are described mathematically by non- 
analytic behaviour of thermodynamic functions, which re- 
flect the drastic changes taking place in the system at a mi- 
croscopic level. Besides materials science, phase transi- 
tions play vital roles in cosmology, particle physics, chem- 
istry, biology, sociology and beyond; The universe began 
in a symmetric manner and went through a series of phase 
transitions through which the particles of matter with which 
we are familiar (electrons, protons, the Higgs boson, etc.) 
materialised. More mundane examples include traffic flow 
(where there is a transition between jammed and free-flow- 
ing states), growth phenomena and wealth accumulation 
(where there may be a transition to a condensation phase, 
for example). While the latter examples refer to non-equil- 
ibrium systems, the emphasis in this article is on a famous 
transition which exists in equilbrium systems. 

The mathematical physics which describes such phe- 
nomena belongs to the realm of equilibrium statistical me- 
chanics, one of the most beautiful, sophisticated and suc- 
cessful theories in physics. Equilibrium statistical physics 
is based on the following premise: the probability that a 



system is in a state S with energy E at a temperature T is 

where (3 = l/fc^T and ks is a universal constant, and 
Zl(J}) is a normalising factor known as the partition func- 
tion, 

Z L ((3)=J2e- 0E(S) - (2) 
s 

Here, the subscript L indicates the linear extent of the sys- 
tem. 

A related fundamental quantity is the free energy, //, (/?), 
given by 

U{p) = ^2 In Z L {0), (3) 

where d is the dimensionality of the system. Phase transi- 
tions can only occur when the system under consideration 
has an infinite number of states in which it can exist - for 
example, in the thermodynamic limit L — > oo. 

In the modern classification scheme, such phase transi- 
tions are categorised as first-, second- (or higher-) order if 
the lowest derivative of the free energy that displays non- 
analytic behaviour is the first, second (or higher) one. Tran- 
sitions of infinite order brake no system symmetries. The 
most famous of these is the Berezinskii-Kosterlitz-Thouless 
(BKT) transition in the two-dimensional XY model |l]|2). 

2. THE TWO-DIMENSIONAL XY MODEL 

The model is defined on a two-dimensional regular lattice, 
whose sites are labeled by the index i, each of which is 
occupied by a spin or rotator si. These two-dimensional 
unit vectors have 0(2) or U(l) symmetry. The energy of a 
given configuration is 

E=-^2 SiSj , (4) 

where the summation runs over nearest neighbouring sites 
or links. This model is used to study systems such as films 
of superfluid helium, superconducting materials, fluctuat- 
ing surfaces, Josephson-junctions as well as certain mag- 
netic, gaseous and liquid-crystal systems. 

The scenario proposed in seminal papers by Berezinskii 
and Kosterlitz and Thouless 1 1 2) is that at a temperature 
above the critical one (T > T c or j3 < (3 C ) positively 
and negatively charged vortices (i.e., vortices and antivor- 
tices) which are present (see Fig. 1) are unbound (dissoci- 
ated from each other) and disorder the system. Below the 
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Figure 1: Examples of charge-1 (left) and charge-2 (right) 
vortices in the XY model on a lattice. Spins rotate through 
2tt and Att, respectively, as contours around the two bold 
points are traversed. 

critical temperature (T < T c or (3 > (3 C ) they are bound 
together and are relevant as dynamical degrees of freedom. 
There, long-range correlations between spins at sites i and 
j (separated by a distance r, say) exist and are described 
by the correlation function whose leading behaviour in the 
thermodynamic infinite-volume limit is 

Goo(r)~r-"P). (5) 

The correlation length (which measures to what extent spins 
at different sites are correlated) diverges and this massless 
low-temperature phase persists, with the system remain- 
ing critical with varying rj((3), up to (3 — f3 c at which 
rj(/3 c ) = rj c . Above this point, correlations decay expo- 
nentially fast with leading behaviour 

Goo(r) ~ e --/e-C*>. (6) 

Here £oo(i) is the correlation length, and t = T/T C — 1 > 
measures the distance from the critical point. 

As this critical point is approached, the leading scaling 
behaviour of the correlation length, the specific heat and the 
susceptibility (which respectively measure the response of 
the system to variations in the temperature and application 
of an external magnetic field) are 

&»(*) ~ e bt ~" , (7) 
Gao(t) ~ C 2 , (8) 
Xo.it) ~ gT", (9) 

in which and & is a non-universal constant. This exponen- 
tial behavour is known as essential scaling, to distinguish 
it from more conventional power-law scaling behaviour (in 
which, for example, £oo(i) ~ t~ u ). In summary, the BKT 
scenario means a transition which (i) is mediated by vortex 
unbinding and (ii) exhibits essential scaling. 

Besides the two-dimensional XY model, transitions of 
the BKT type exist in certain models with long-range in- 
teractions 1 3 1, antiferromagnetic models |4|, the ice-type F 
model |5| and in string theory |6| amongst others. Thus 
a thorough and quantitative understanding of the paradig- 
matic XY model is crucial to a wide breadth of theoretical 
physics. 

For many years Monte Carlo and high-temperature anal- 
yses of the XY model sought to verify the analytical BKT 
renormalisation-group (RG) prediction that v = 1/2 and 



rjc = 1/4 and to determine the value of (3 C . Typically 
(3 C was determined by firstly fixing v = 1/2. Subsequent 
measurements of i] c yielded a value incompatible with the 
the BKT prediction. Because of the elusiveness of its non- 
perturbative corroboration, the essential nature of the tran- 
sition was questioned |7|. See Table 1 of |8| for an exten- 
sive overview of the status of the model up to that point. 

3. LOGARITHMIC CORRECTIONS 

The fundamental realization of £0 was that the thermal 
scaling forms ^ and l|9} (and similar formulae for related 
thermodynamic functions) are inconsistent as they stand. 
Instead, they must be modified to include logarithmic cor- 
rections: CooO) ~ £oo(t)~ 2 {\n^ 00 (t)) q with q = 6 and 

Xoo(t) ~ M*) 2 -"' (ln£ooWr 2r , (10) 

where RG indications implicit in 1 2 1 are that r = — 1/16 = 
-0.0625. 

Eq.(I10> is an analytic prediction, which, since it is based 
on perturbation theory, requires confirmation through non- 
perturbative approaches. However, infinite lattices are not 
attainable using finite numerical resources. Instead, one 
simulates finite systems, where, at criticality (t = 0), the 
lattice size L plays the role of ^ in this model. The finite- 
size scaling (FSS) prediction for the susceptibility is then 

XL(0)~L 2 -^(lnL)- 2r . (11) 

To verify the BKT scaling scenario then, the thermal for- 
mula (II Oi and/or the FSS formula il 1> needs to be con- 
firmed numerically. Sophisticated FSS techniques (involv- 
ing partition function zeros) were used in JH] |9| to resolve, 
for the first time, the hitherto conflicting results for v, j3 c 
and rye. However, the analysis resulted in an estimate of 
—0.02(1) for r, a value in conflict with the RG prediction 
of r = -1/16 = -0.0625 from 0. 

Thus, in recent years, the focus of numerical studies of 
the XY model shifted to the determination of the logarithm 
exponent r. Indeed, the FSS analyses of | lOl ll lllT2l . using 
(lilt , yielded values compatible with that of 1 8 9 1 but in- 
compatible with |2| (see Table 1). Nonetheless, it was clear 
that taking the logarithmic corrections into account leads 
to the resolution of the v-(3 c -rj c controversy. The most pre- 
cise estimate for the critical temperature in the literature for 
the XY model with the standard action of (0} is contained 
in HI 31 and is (3 C = 1.1199(1). This value was obtained 
by mapping the XY model onto the exactly solvable body 
centered solid-on-solid model, and thereby circumventing 
the issue of logarithmic corrections. As demonstrated in 
Table 1, recent analyses which have inluded the logarith- 
mic corrections have resulted in estimates for f3 c compati- 
ble with this value (while keeping r\ = 1/4). In (Tl | and 
111 21 respectively, phase transitions in a lattice grain bound- 
ary model and a lattice gauge theory are studied. Because 
these are in the same universality class as the XY model in 
two dimensions, they have the same scaling behaviour (lilt . 
However, the critical temperatures in these models bear no 
relationship to that of the XY counterpart. Although 1 10 1 
contains a study of the XY model, it employs the Villain 



Table 1: Estimates for j3 c and r for the XY model from 
a selection of recent papers with indications of the method 
used to obtain them (RG related, FSS or thermal scaling). 



Authors 


Year 


Method 


0c 


r 


Kosteiiitz, 


1973 


RG 




-0.0625 


Thouless (2) 










Irving, 


1997 


FSS 


1.113(6) 


-0.02(1) 


Kenna 










Patrascioiu, 


1996 


thermal 




0.077(46) 


Seiler fl4| 










Campostrini 


1996 


thermal 


1.1158(6) 


0.042(5) 


et al. (15) 






1.120(4) 


0.05(2) 


Janke tlOl 


1997 


FSS 
thermal 




-0.027(1) 
0.0560(17) 


Hasenbusch, 


1997 


RG 


1.1199(1) 




Pinn Il3l 










Jaster, 


1999 


FSS 




-0.0233(10) 


Hahn fTTI 




thermal 




0.056(9) 
0.070(5) 


Dukovski 


2002 


FSS 


1.120(1) 




et al. fl6l 










Chandrasekharan, 


2003 


FSS 




-0.035(10) 


Strouthas fj2) 










Hasenbusch 1171 


2005 


FSS 




-0.056(7) 



formulation, which is different to (TJJ, and leads to a differ- 
ent value for the nonuniversal quantity (3 C . 

It was suggested in 1 8 10 1 that just as taking the log- 
arithmic corrections into account leads to the resolution of 
the leading scaling behaviour, in the same spirit it is con- 
ceivable that numerical measurements for r may become 
compatible with the RG prediction if sub-leading correc- 
tions are taken into account. In this case, (lilt is more fully 
expressed as 

Xi (0)^L 2 --(lnLr 2 "{l + o(^)} . (12) 

This was tested in 1 8 10 1 to determine if the discrepancies 
between the numerical and theoretical estimates for r can 
be ascribed to sub-leading corrections. However, the lattice 
sizes available (up to L = 256) were too small to resolve 
the issue. 

Recently this issue was again addressed in a very high 
precision numerical simulation (using lattices as big as L = 
2048) in ATI . Using the ansatz 

X L (0)~L 2 -^(C + lnL)- 2r , (13) 

and fitting to r gives r — —0.056(7) if large enough lat- 
tices are used. This value is compatible with the analytic 
prediction, r = —0.0625. 

Notwithstanding this result, which is FSS based, it is 
rather surprising that all of the analyses of the thermal scal- 
ing formula i ll Ob in Table 1 yield positive values of r far 
from the RG prediction that r = —0.0625. It appears, 
therefore, that the FSS approach is more powerful than that 
based on thermal scaling, although more extensive analy- 
ses would be required to test how this approach compares 
with that of fTTI. 



While large scale simulations were required to resolve 
these puzzles, a different technique, not requiring extensive 
simulations was used in 1 1 8 1 . A conformal mapping was 
used to switch from a confined L x L lattice to the semi- 
infinite half plane. In this way, ?y(/3) could be deduced from 
the congelation function at any value of (3 > (3 C . In par- 
ticular, 7/ c = 1/4 is accurately recovered at the transition 
temperature and clear evidence for existence of logarithmic 
conections in the conelation function G(r) is presented. 

4. VORTEX UNBINDING 

The vortex-binding scenario is crucial to the BKT phase 
transition in the two-dimensional XY model. Because the 
energy of a single vortex increases with the system size 
as In L, at low temperature they can only occur in vortex- 
antivortex pairs. Mutual cancelation of their individual or- 
dering effects means that such a pair can only affect nearby 
spins and cannot significantly disorder the whole system. 
Topological long-range order exists in the system at low 
temperature. However at high temperature, the number 
of vortices proliferates and the distance between erstwhile 
partners becomes so large that they are effectively free and 
render the system disordered. 

It was for a long time believed that altering the ener- 
getics of the XY model to disable the vortex-binding sce- 
nario may lead to a transition different to the Berezinskii- 
Kosterlitz-Thouless one |19|. The step model is obtained 
from the XY model by replacing the Hamiltonian © by 

£ = -^sgn(s i s i ) . (14) 

The energy associated with a single vortex for this system 
is expected to be independent of the lattice size. There- 
fore, on the basis of this argument, vortices could exist at 
all temperatures, disordering the system. I.e., there was ex- 
pected to be no vortex-driven phase transition in the step 
model - if there is a phase transition, it was expected not 
to be of the BKT type. Indeed, early studies supported this 
assertion 1 19 1. 

However, in 1 8 , 9 1, very strong numerical evidence was 
presented that (i) there is a phase transition in the step model 
and (ii) it is of the BKT type (with even the conections to 
scaling being the same as those for the XY model). Given 
the very different vortex energetics for the two models, this 
came as a surprise. 

The issue was further addressed in 1 20 1, where evidence 
for the existence of a BKT transition in the step model was 
again proffered. The approach of |20| focused on numer- 
ical analyses of the helicity modulus, which experiences 
a jump at the transition. A similar approach to the XY 
model is contained in 1211 . The main idea of 1201 . which 
explains the occurence of the BKT transition in the step 
model, is that while the energy associated with a single 
fixed vortex in the system remains finite, the free energy 
grows as In L. This fact inhibits proliferation of free vor- 
tices in the low-temperature phase. It further implies that 
the harmonic properties of the interaction @ do not form 
a necessary condition for a BKT transition. Consequently, 
and as pointed out in |20|, the BKT phase transition may 



be an even more general phenomenon than hitherto recog- 
nised. 

5. FUTURE DIRECTIONS 

Asymptotic freedom of d — 2 O(N) models: It is gen- 
erally believed that there are differences of a fundamental 
nature between abelian and nonabelian models. The XY 
model has an 0(2) symmetry group and is abelian, while 
all 0(N) models with N > 2 are nonabelian. The Mermin- 
Wagner theorem states that a continuous symmetry of the 
O(N) type cannot be broken in two dimensions |22|. Thus 
there cannot be a transition to a phase with long-range order 
in either of the N = 2 or N > 2 scenarios there. However, 
in a two-dimensional theory, topological defects of dimen- 
sion rn can exist if the (1— m) th homotopy group, 7Ti_ m , of 
the order parameter space is non-trivial. For O(N) models, 
this space is the hypersphere S 1 ^^ 1 . The only non-trivial 
group is tti(S 1 ) which is isomorphic to the set of integers 
under addition. This is the condition that gives rise to point 
defects (vortices) with integer charge in the N = 2 case 
(the XY model). The binding of these vortices at low tem- 
perature is the mechanism giving rise to the BKT phase 
transition 1 l]|2). 

For N > 2, conditions are not supportive of the ex- 
istence of topological defects of this type and the widely 
held belief is that there is no phase transition, there being 
no distinct low-temperature phase. Perturbation theory pre- 
dicts that the N > 2 models are asymptotically free. There 
is, however, no rigorous proof to this effect. This belief 
has been questioned and 1 14 23 1 have given numerical ev- 
idence for the existence of phase transitions of the BKT 
type in these models as well as heuristic explanations of 
why such transition could occur and a rigorous proof that 
this would be incompatible with asymptotic freedom. 

Perturbative and Monte Carlo calculations for the N — 
3 and N = 8 models have been performed in II24I . which 
do not support the existence of such a BKT-like phase tran- 
sition there and instead are in agreement with perturbation 
theory and the asymptotic freedom scenario. Nonetheless, 
the controversy has not entirely gone away 1251 1261 1271 , 
and one may argue that inclusion of logarithmic considera- 
tions could help for a precise unambiguous resolution. 

The diluted XY model: An interesting current topic of 
research has been the question of the role of impurities in 
the XY model |28|. The presence of impurities brings the 
model closer to real systems, where such physical defects 
are present. Impurities are modeled by randomly diluting 
the number of sites (or bonds) on the lattice. Clearly, if the 
dilution is so strong as to inhibit the percolation of spin- 
spin interactions across the lattice (such that it is effectively 
broken into finite disconnected sets) no phase transition can 
occur for any model. Thus moderate dilution generally is 
expected to decrease the location of the transition temper- 
ature. However, the special additional feature of the XY 
model is the presence of vortices and the fact that they drive 
the transition. Vortices are attracted to and, to some extent, 
anchored by impurities and the vortex energy is reduced at 
such a vacancy. Therefore with increasing dilution, more 
vortices can be formed and the amount of disorder in the 



system is increased. This effect may enhance the lowering 
of the critical temperature to such an extent that it vanishes 
before the percolation threshold is reached. 

This is the issue addressed in 1 28 29 30 31]. The per- 
colation threshold occurs when the density of site vacancies 
is p w 0.41. The critical vacancy density is identified by 
the vanishing of the critical temperature, which in turn is 
identified as being the location at which r/(f3) — 1/4. In 
1 29 1 the critical temperature was reported to vanish above 
the percolation threshold at vacancy density p w 0.3. How- 
ever, in 1301 it was suggested that the critical density is, in 
fact, closer to the percolation threshold. Support for the lat- 
ter result recently appeared in 1311 . If this is true, it means 
that the vortices do not, in fact, strongly enhance the low- 
ering of the critical temperature. 

From the collective experience with the pure XY model 
as reported above, it is clear that a more precise identifica- 
tion of the critical temperature through r)(0) = 1/4 would 
require taking account of the logarithmic corrections, al- 
though ignoring them may suffice as first approximation. 
Indeed, ignoring these corrections in the pure model also 
leads to an estimate for the critical temperature which is 
higher than accurate values 1 3 1 1. 

Besides the value of the critical temperature in a diluted 
model, one is also interested in the scaling behaviour of the 
thermodynamic functions at the phase transition there. The 
Harris criterion predicts that disorder does not change the 
leading scaling behaviour of a model if the critical expo- 
nent a associated with the specific heat of the pure model 
is negative |32|. This is the case for the XY model in 
two dimensions. However, it is unclear what effect dilution 
can have on the quantitative nature of the exponents of the 
logarithmic corrections in such a case. This would be an 
interesting avenue for future research, and the XY model 
(which has negative a) offers an ideal platform upon which 
to base such pursuits. 

Other models with logarithmic corrections: Logarithmic 
corrections to scaling also exist in other important mod- 
els. While their existence has been unambiguously ver- 
ified in d = 4 O(N) models J33] |34|, this is not so in 
most cases |35|. In particular, Kim and Landau applied 
FSS techniques to the four-state d = 2 Potts model |36|. 
Again, FSS behaviour could only be described if the loga- 
rithmic corrections are included. However, here inclusion 
of the leading multiplicative logarithmic corrections is in- 
sufficient and sub-leading additive corrections are required. 
The puzzle of why unusually large numbers of correction 
terms are necessary (despite the availability of the exact 
value of (3 C in this model) could, perhaps, now be resolved 
by an analysis on the scale of [ 17 1 and this would be another 
interesting avenue to pursue. Similar problems concerning 
logarithmic corrections and their detection exist in other 
two-dimensional models such as the diluted Ising model 
11371 . Theoretical progress on the general issue of logarith- 
mic corrections will be reported elsewhere |35 1. 

6. CONCLUSIONS 

The issue of topologically driven phase transitions charac- 
terized by the two-dimensional XY model has been revis- 



ited and a timely review of the status of scaling at the fa- 
mous Berezinskii-Kosterlitz-Thouless transition given. Af- 
ter two decades of work, the perturbative renormalization 
group precictions of [ 1 2 1 for the leading scaling behaviour 
were confirmed in 1 8 9 1 and there is now little or no doubt 
about the correctness of the analytic predictions, v = 1/2 
and T] c — 1/4. The recent controversy over the value of 
the leading logarithmic correction exponent r (summarized 
in Table 1) has been re-examined and claims as to its res- 
olution (at least in the context of finite-size scaling) and 
state-of-the-art calculations summarized 1 17 1. Besides the 
XY model, such multiplicative logarithmic corrections are 
manifest in a variety of other models, and their resolution 
in these contexts is now at the forefront current modern nu- 
merical investigations in statistical and lattice physics. A 
summary of some recent work on such models, as well as 
likely future directions has been given. 

Finally, recent work 1 20 1 confirming the vortex-binding 
scenario as the phase transition mechanism in the XY model 
has also been reviewed to complete a full account of the 
current status of one of the most remarkable and beautiful 
models in theoretical physics. 
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